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Abstract 

Thermal effects are incorporated into previously developed discrete layer mechanics for 
piezoelectric composite beam structures. The updated mechanics explicitly account for 
the complete coupled thermoelectromechamcal response of smart composite beams. 

This unified representation leads to an inherent capability to model both the sensory and 
actuator responses of piezoelectric composite beams in a thermal environment. Finite 
element equations are developed and numerical results are presented to demonstrate the 
capability of the current formulation to represent the behavior of both sensory and active 
smart structures under thermal loadings. 


Introduction 


The development of intelligent and smart composite materials with piezoelectric compo- 
nents offers great potential for use in advanced aerospace structural applications. This 
potential arises from two basic characteristics of piezoelectric materials, called the direct 
and converse piezoelectric effects, which allow for their use as sensors and actuators in 
composite structures. In the direct piezoelectric effect, application of a mechanical force 
on the piezoelectric material leads to the induction of an electrical charge or voltage, 
while in the converse piezoelectric effect the application of a charge or voltage produces 
a mechanical force or strain in the piezoelectric material. By taking advantage of the 
direct and converse piezoelectric effect of embedded piezoceramic or piezopolymer 
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devices, these novel materials can combine the superior mechanical and thermal proper- 
ties found in conventional composites along with the inherent self-monitoring and 
adaptive capabilities of piezoelectric materials. An area where such materials may 
provide dramatic advantages is in the development of smart composite structures with 
the capability to sense thermally induced distortions and to actively compensate for 
adverse thermomechanical conditions. Typical applications of such structures are 
envisioned in the thermal distortion management of propulsion components and/or space 
structures. However, before they can be utilized in these applications the performance of 
smart piezoelectric structures in thermal environments must be quantified. Consequently, 
this paper will present the development of comprehensive mechanics for the analysis of 
such smart thermal composite structures. 

Extensive development of analytical methods for modeling the behavior of piezoelectric 
structures has been reported in the literature. The early works focused mainly on the 
development of simplified models for the vibration control of beams [Bailey and Hub- 
bard (1985), Crawley and deLuis (1987)]. This was followed by the extension of classical 
lamina tion theory to incorporate piezoelectric effects [Tzou and Gadre (1989), Lee and 
co-workers (1991), Wang and Rogers (1991)]. The development of this piezoelectric 
laminate theory lead in turn to the formation of finite element methods for the analysis 
of smart structures using plate and shell elements [Allik and Hughes (1970), Naillon and 
co-workers (1983), Tzou and Tseng (1990), Lammering (1991), Ha and co-workers 
(1992)]. However, since classical lamination theory has been shown to have certain 
limitations in the analysis of thick or heterogenous (through-the-thickness) laminates 
[Reddy (1987)], a discrete layer theory for active beams based on an induced strain 
approach was implemented by Robbins and Reddy (1991) into a finite element formula- 
tion. This discrete layer theory was in turn extended to incorporate the coupled equa- 
tions of piezoelectricity and developed into a multi-field formulation for beam and plate 
finite element equations by Heyliger, Ramirez & Saravanos (1993) and Saravanos & 
Heyliger (1994). 

All of the previously described methods neglect the implications of thermal effects on 
both the active and sensory response of smart structures. Thermal effects become 
important when the piezoelectric structure has to operate in either extremely hot or cold 
temperature environments. These extreme conditions may severely affect the response 
of smart structures in three distinct ways: (1) induction of thermal stresses due to 
differences in the coefficients of thermal expansion between the various composite plies 
and piezoelectric layers, (2) pyroelectric phenomena (i.e. changes in the electrical 
displacements arising from the coupled electrical and thermal behavior of the piezoelec- 
tric material), and (3) temperature dependence of the elastic, piezoelectric, and dielectric 
properties. Only a limited amount of work has been reported concerning this topic. 
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Mindlin (1974) formulated a set of two-dimensional thermopiezoelectric equations for 
plates. Tauchert (1992) extended classical lamination theory to account for the coupled 
effects of thermopiezoelectricity. Rao and Sunar (1993), Tzou and Howard (1994), and 
Tzou and Ye (1994) developed finite element formulations for plates and shells (based 
on classical lamination theory) to account for thermopiezoelectricity. 

Consequently, the purpose of this paper is to extend the previously developed discrete 
layer formulation of Saravanos and Heyliger (1994) to account for the coupled mechani- 
cal, electrical, and thermal response in modern smart composite beams. The mechanics 
account for thermal effects which may arise in the elastic and piezoelectric media at the 
material level through the constitutive equations. The displacements, electric potentials, 
and temperatures are introduced as state variables, allowing them to be modelled as 
variable fields through the laminate thickness. This unified representation leads to an 
inherent capability to model both the active compensation of thermal distortions in smart 
structures and the resultant sensory voltage when thermal loads are applied. The 
corresponding finite element formulation is developed, and numerical results demonstrate 
the ability to model both the active and sensory modes of composite beams with 
heterogeneous plies with an attached piezoelectric layer under thermal loadings. 


Governing Material Equations 

This section outlines the foundations and steps required for developing the governing 
equations of composite beams with embedded piezoelectric layers. The mechanical 
response of the piezoelectric material can be represented by the equation of motion 

o u. - O- ■ + f- 

r i yj Ji 


where p, u i , u ; , a,, f. are the density, acceleration, displacement, stress, and body 
force per unit volume, respectively; with i, j ranging from 1 to 3. 

The electrical response of the piezoelectric material is described by the electrostatic (or 
Maxwell’s) equation 

D u = 0 < 2 

where Dj is the electric displacement. 

The small deformation strain-displacement relations are 
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(3) 


5 * 2 U 


V + U u 


where is the strain tensor. The electric field vector E 4 is related to the electric 
potential <p by 

Et = - 


(4) 


In the lin ear theory of piezoelectricity, the equations of elasticity are coupled to the 
equations of electrostatics. The constitutive equations for a thermopiezoelectric material 
[Nye (1964)] employing standard contracted notation are 



S. - s.Y(n o„ * <£( T) E. ♦ o Y(T) e 

(5) 


D m - <£(D E t * P : r (T) e 

(O 

or in semi-inverted form 

s t - <£(d e, - xYm e 

(7) 


D m *'I.(T) s, * 4 r (D e, ♦ pO; t (T) e 

(8) 

with, 

*L.<t) = cff(T) dl.cn 
xf r (D - c.y<d «.fcn 

0 = T - T 0 

(9) 


<g(T) - ei r (D - 4(7) 4(D 
p s . T cn - P ;- T cn - <&cn 



where D,,, is the electric displacement component; and s^ are the elastic stiffness and 
compliance tensors; d^ and e^, are the different forms of the piezoelectric tensor; e mk is 
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the electric permittivity tensor; a a and k a are the different forms of the coefficient of 
thermal expansion; p m is the pyroelectric constant; 0 is the temperature difference from 
the stress free reference temperature T c to the current temperature T; superscripts E, o, 
S, and T, represent constant voltage, constant stress, constant strain, and constant 
temperature conditions, respectively; a, p = 1, ..., 6; and k, m = 1, 2, 3. 

Through use of the divergence theorem and neglecting body forces, Eqs. (1) and (2) can 
be expressed in an equivalent variational form as 

f v P & i bu i dv + f v °i> bS v dV ' f v , bEi dV (10) 

= dA + q 6<J> dA 

where ~t i are the surface tractions applied on the surface A; q is the electrical charge 

applied on the surface Ap of the piezoelectric material; V represents the whole volume 
including both composite and piezoelectric materials; and V p represents the volume of 
the piezoelectric material only. 


Discrete leaver Formulation 


A discrete layer lamina te theory for thermal piezoelectric composite beams is described 
in this section. For beams, a simplified displacement and electric potential can be 
obtained since only axial variations of the state variables are assumed 

u = u(x,z ) , V = 0 , w = w(x,z) , 4> = 4>(*>Z) (U 


where x is the in-plane axis of the beam and z is the through-the-thickness axis. 

The discrete layer theory implements the following piecewise continuous approximations 
through the laminate thickness for the state variables 


u(x,z,t ) = £ U J (x,t ) V(z) 
j - 1 


( 12 ) 


w(x y z,t) = £ W*{x,t) V w (z) 


(13) 
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(14) 


J» 


♦<w) = E ♦W) Vu) 


e(z,o = £ 0^(0 V(z) 


/-i 


(15) 


where n is the number of interpolation functions For the case of constant through- 
the-thickness displacement [ w(x,z,t) = w°(x,t) ], is equal to unity. Currently, i|r J is 
represented using linear Lagrangian interpolation functions. The simplified displace- 
ments and the assumption of constant through-the-thickness displacement leads to the 
simplified constitutive equations derived from Eqs. (5)-(8) shown in Appendix A. 

By combining Eqs. (3), (7), (8), (10), (12)-(15) and integrating through-the-thickness, the 
following variational form can be obtained for the case of a constant through-the- 
thickness displacement beam 


8k 1 dx + J P 33 w 6w dx 


k - 1 m-1 


n n 


EE / < D n 


* du - as** „ D * u . tul ) dx . 


k* 1 *»*1 


cbc dx 


fa j 55 dx dx J 35 dx dx 


X 


Jr -1 n -1 


dx dx 


( 16 ) 


EE / e£ <t>- a* 1 dx - E / </i* - -«*»+*) 

Jt-l «-l *-l * 


= E t/f «***«* «♦*].♦[/, »*>], 


k-l 
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The discrete layer laminate matrices for the density, stiffness, piezoelectric, electric 
permittivitty, and external forces of a composite beam have been presented previously in 
a general form by Saravanos and Heyliger (1994) and are included for the specific case 
of a constant through-the-thickness displacement beam in Appendix B. The newly 
derived resultant laminate thermal force vector is 

/£* = £ / b *1 Vc*) ♦*<*) * (17) 

1-1 Jz 

and the resultant thermal electric displacement is 

q& = £ / b P 3 & V(z) d* (18) 

where b is the width of the beam and L is the number of plies. 


Finite Element Formulation 


The finite element formulation for a composite piezoelectric beam is obtained by 
incorporating additional local in-plane approximations to the state variables 


m n 


u(x,z,t ) = £ £ UHt) R\x) V(z) 
1*1 7*1 


m it 


w(x,z,f) = £ £ w*(f) *'(*) V(z) 

i*l 7*1 


= £ £ V'(0 *'(*) V(z) 

i-i ;-i 


( 19 ) 

( 20 ) 
( 21 ) 


where m is the n um ber of in-plane shape functions R. Currently, R is represented using 
linear Lagrangian interpolation functions and a selectively reduced integration scheme is 
implemented for the second stiffness term (containing D 55 ) of Eq. (16) in order to 
eliminate overstiffening at low thicknesses. 


By combining Eqs. (3), (7), (8), (10), (19)-(21), the following finite element matrix 
formulation is obtained for the case of a constant through-the-thickness displacement beam 
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[M n ] 0 0 


{#} 


[K lt ] [K 13 ] [K u ] 


{*/} 

0 [Hf 33 ] 0 

1 

{W} 

► + 

[*,,] [**] 0 

< 

{W } > 

0 0 0 




[*4il 0 [K u ] 


.{*} . 


W 



{ f 3 ) 

► + « 

0 

. . 


{<?*> 


( 22 ) 


where the elements of the sub matrices above are calculated from the generalized 
discrete layer laminate matrices defined previously in Appendix B and in Eqs. (17-18). 
The finite element matrices for the mass, stiffness, and external force have been present- 
ed in a general form by Saravanos and Heyliger (1994) and are included in Appendix C 
for the specific case of a constant through-the-thickness displacement beam. The newly 
derived finite element thermal force vector is 

FU - lu ^ * (23) 

and the thermal electric displacement is 

Qi, - - £ 4* R'w dx (24) 


The coupled finite element formulation can also be expressed in a compact form as 


[MJ 0 

4 

f«*»L 

[*J [KJ 

4 

(i/n . f(Fi 



o o 


[ml 

[**] [K„] 


lm I * tail 




The electric potential vector in Eq. (25) can be further partitioned into active and 
sensory components such that {0} = {0 a ; 0 s } in order to obtain 
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(26) 


[MJ 0 
0 0 


{U} 


I*J [*S3 

[*£] [*53 


{*‘1 


{F} + {F m } - [K%] {**} 

[Q s ] * {Q*} - [^3 . 


where the superscripts s and a indicate the partitioned submatrices in either sensory or 
active configuration respectively. This form has the advantage of positioning the 
unknown variables (displacements and sensory electric potentials) in the left-hand terms, 
while the known quantities (mechanical loads, thermal loads, electric charges, and active 
voltages) are included in the right-hand terms. 

Eq. (26) can then be uncoupled into the following independent equations for the 
structural displacements 

[M„] tui * ( [XJ - wjr 1 [*£] ) (W) * (FI * (27) 

(f.i * ( [*£] [*£]-' uQ - nQ > «**> 

and the electric potentials at the sensors 

{9 s } = -[JQ- 1 ( [*£] {U} + [K£] {*«} - {<?*} - { Q *} ) (28) 


A pplications 


Materials and Assumptions 

Results from representative problems are presented in this section to evaluate the 
discrete layer formulation and to demonstrate the performance of simple sensory/active 
structures in thermal environments. Two different composite cantilever beam configu- 
rations are considered as shown in Figures 1(a) and 1(b). For both cases, the beam is 
25.4 cm (10 in) long and consists of various orientations of graphite/epoxy plies attached 
to piezoelectric layers, where each graphite/epoxy ply and piezoelectric layer has a 
thickness of 0.0127 cm (0.005 in). A total of 30 equally spaced linear beam elements are 
modeled along the length of the beam, while one discrete layer is used through-the- 
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thickness for each piezoelectric and graphite/epoxy layer as shown in Figure 2(a). The 
room temperature material properties used for the graphite/epoxy and piezoelectric 
layers are listed in Table I. The piezoelectric material used for this study is the commer- 
cially available piezoceramic APC 840 (American Piezo Ceramics Inc.), while representa- 
tive material properties of a graphite/epoxy composite were selected. The current study 
uses the constant through-the-thickness displacement beam theory and neglects the 
temperature dependence of material properties. 

[fl q/p] Beam Under Uniform Thermal Load 

Thermal Displacements. The first beam considered is a [Og/p] configuration, where p 
represents the piezoelectric layer, under a uniform thermal load of 100 °C. The presence 
of the piezoelectric layer results in an asymmetric laminate configuration, which induces 
thermal distortions under the thermal load. Figures 3 and 4 illustrate the induced axial 
and transverse displacements produced when the beam is operating in closed circuit 
conditions. The closed circuit condition is obtained by grounding (0 Volts) both the 
lower and upper surfaces of the piezoelectric layer. Grounding both surfaces effectively 
eliminates the induced piezoelectric strains and produces a virtually conventional 
composite beam. Thus, the results for this grounded loop condition can also be validated 
using a conventional finite element structural analysis. This conventional finite element 
analysis was performed using the MARC analysis code (MARC Analysis Research 
Corporation) with nine layers of 15 eight-noded quadratic plane stress quadrilateral 
elements as shown in Figure 2(b). Both methods result in almost identical axial and 
transverse thermal displacements, as shown in Figures 3 and 4, lending credence to the 
accuracy of the current formulation. 

The resultant thermal deflections of the [Og/p] beam under the open circuit or sensory 
conditions (only the lower surface of the piezoelectric layer is grounded) are shown in 
Figures 5 and 6. Since only the lower surface is grounded, the effects of the piezo- 
electrically induced strains are no longer neutralized and leads to the formation of a 
sensory electric potential on the top surface as a result of the induced thermal stresses. 
Figure 6 shows a significant difference (~15 %) in the transverse displacement between 
the closed and open circuit conditions. This demonstrates the importance of including 
the thermo-electrical coupling to model the behavior of thermal piezoelectric structures 
and illustrates the ability of the mechanics to capture this phenomena. 

Sensory Voltage. The open loop case also demonstrates an important potential applica- 
tion of sensory structures in the monitoring of thermal distortions. In these applications, 
the displacement fields which arise under thermal loadings may be sensed through 
correlation with the measured voltages at the sensors. This capability is demonstrated in 
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Figure 7, which shows the measured electric potential (from both the thermal strain and 
pyroelectric effects) on the top surface corresponding to the axial and transverse 
displacements shown in Figures 5 and 6. The thermal strain sensory voltage accounts for 
the effects arising from the thermally induced deformations, while the pyroelectric voltage 
quantifies the effects that occur due to the variation in temperatures. This application 
demonstrates the capability to monitor the thermal displacements via the different 
sensory voltages, which will provide the inferred essential feedback required for thermal 
distortion management of smart composite structures. 

Thermal Distortion Management. The capability to actively compensate for thermal 
distortions is demonstrated in Figure 8 for the [Og/p] beam. Through application of a 
voltage differential between the upper and lower surfaces of the piezoelectric layer, the 
converse piezoelectric effect can be utilized to compensate for the thermal displacements. 
As shown in the figure, the thermal displacement behavior of the beam can be complete- 
ly altered from the open loop state depending on the value of the voltage differential. 

The thermal distortion can be eliminated by applying a voltage differential of 300 V, 
while a completely opposite curvature of the beam can be obtained by applying a voltage 
differential of 500 V. These results indicate the potential of piezoelectric structures to 
suppress thermally induced bending. 

[p/0/90/±45] Ream Under Thermal Gradient 

The present e xam ple considers multiple piezoelectric layers (to simultaneously take 
advantage of both the sensory and active behavior of piezoelectric materials), a more 
general lamination [p/0/90/±45] s , and a through-the-thickness linear thermal gradient 
(120 °C on the top surface and 20 °C on the bottom surface) as shown in Figure 1(b). 
Even though the beam has a symmetric laminate configuration including the piezoelectric 
layers, the applied thermal gradient produces a bending deflection in the beam. 

Thermal Displacements. The induced axial and transverse displacements produced in the 
beam when both piezoelectric layers are operating in closed circuit conditions is shown in 
Figures 9 and 10. Also shown are the results from a conventional finite element analysis 
performed using the MARC analysis code (MARC Analysis Research Corporation) with 
ten layers of 15 eight-noded quadratic plane stress quadrilateral elements. Both methods 
lead to identical predictions, lending credence to the validity to the current formulation. 


Combined Sensory and Actuator Application. In this application, the upper piezoelectric 
layer is assumed to act as the sensor to monitor the induced electric potential, while the 
lower piezoelectric layer is used to actively compensate thermal deflections. Figures 11 
and 12 show the axial and transverse deflections of the beam under applied voltages of 0 
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V, 100 V, and 200 V on the active (lower) piezoelectric layer. As shown in the figures, 
increasing the applied voltage reduces the thermally induced deflections from the 
reference state (0 V). The corresponding sensory electric potentials (from both the 
thermal strain and pyroelectric effects) measured on the upper piezoelectric layer 
(sensor) for the different active voltages are shown in Figure 13. In actual applications, a 
control loop would be established between the sensory and active layer, to monitor the 
sensory voltages and then to generate the appropriate active voltage which will minimi ze 
the thermal distortions. This capability is planned to be investigated in the near future. 

Convergence Studies 

The accuracy of the beam element developed in this study is verified for both beam 
configurations examined. The open circuit condition case is used for the [0g/p] beam, 
while the 100 V active voltage case is used for the [p/0/90/±45] 8 beam. The results are 
shown in Table II and indicates the relative insensitivity of this element to the varying 
aspect ratios for the configurations examined. 

Summary 

The development of discrete layer mechanics to model the coupled thermoelectro- 
mechanical response of smart thermal composite beams with embedded piezoelectric 
layers was described. The corresponding finite element formulation was presented and a 
linear layerwise beam element developed. The capability of the updated mechanics to 
capture the response of piezoelectric composite structures under a thermal environment 
was demonstrated through numerical examples performed on a graphite/epoxy beam with 
attached piezoelectric layers. For the numerical study, a constant through-the-thickness 
displacement beam theory was used and the temperature dependence of material 
properties neglected. Predicted results for the beam in closed circuit conditions were 
validated using a conventional finite element analysis. 

Applications of the active and sensory modes operating independently were demonstrat- 
ed on a [0g/p] beam under a uniform thermal load of 100* C, while the combined active 
and sensory operations of a [p/0/90/±45] s beam under a thermal gradient was illustrated. 
The numerical studies indicate the significance of thermal effects for smart piezoelectric 
structures and illustrate the strong coupling between thermal and piezoelectric effects 
under open circuit conditions. Overall, the case studies have demonstrated the capabili- 
ties of the current mechanics to represent the complicated thermo-mechanical response 
of thermal smart structures. 
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Nomenclature: 


[AH, [BH, PH 


iq 

[d] 

{D} 


[e] 

{E} 

[E nu 

{f} 


generalized laminate matrices 
elastic stiffness 
piezoelectric tensor 
electric displacement vector 
piezoelectric tensor 
electric field 

generalized laminate piezoelectric matrices 
generalized laminate thermal force vector 




{F} 

[G nm ] 

[K] 

[M] 

P 

[pmn] 

{q> 

q 

{Q} 


body force per unit volume 
thermal force vector 

generalized laminate electric permittivity matrices 

stiffness matrix 

mass matrix 

pyroelectric constant 

generalized laminate mass matrices 

generalized laminate thermal electric displacements 

surface electric charge 

thermal electric displacements 
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R 

[s] 

{S} 

t 

{ t } 

T 

u, v, w 
x, y, z 

XpX^Xa 

a 

S 

M 

X 

♦ 

<t> 

p 

{a} 

0 


in-plane interpolation functions 
elastic compliance 
engineering strain 
time 

surface tractions 
temperature 

displacements along x, y, z axes, respectively 
structural axes 
material axes 

strain coefficient of thermal expansion 

variational operator 

electric permittivity tensor 

stress coefficient of thermal expansion 

through-the-thickness interpolation function 

electric potential 

density 

mechanical stress 

temperature difference from reference temperature 
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Superscripts: 


E 

S 

T 

a 

a 

s 


constant electric field conditions 
constant strain conditions 
constant temperature conditions 
constant stress conditions 
active 
sensory 


Subscripts: 

1 

o 

UU, U0, 00 


layer or ply 
reference conditions 

structural, piezoelectric, and dielectric, respectively 
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Table I. - Room Temperature Material Properties of 
Piezoelectric Composite Beam 


Property 

Graphite/Epoxy 

Piezoceramic 
APC 840 (APC, Inc.) 

Eternity, p, g/cm 3 (lb/in 3 ) 

2.1 (0.076) 

7.60 (0.27) 

Elastic Modulus, E u , GPa (Msi) 

39 (5.7) 

68 (10.0) 

Elastic Modulus, 'E 22 , GPa (Msi) 

8.6 (1.24) 

— 

Shear Modulus, G 12 , GPa (Msi) 

3.8 (0.54) 

26.2 (3.84) 

Major Poisson’s Ratio, v 12 

0.28 

0.30 

Minor Poisson’s Ratio, v 21 

0.06 

— 

Thermal Expansion Coefficient, 
a n , lO^/'C (10“ 6 /°F) 

7.0 (3.9) 

3.8 (2.1) 

Thermal Expansion Coefficient, 
a 22 , lO^-C (10' 6 /°F) 

21.0 (11.7) 

— 

Piezoelectric Charge Constant, 
djj, 10' 12 mA^ 

— 

-125. 

Electric Permittivity, e 33 , 
10- 9 NA' 2 

— 

11.06 

Pyroelectric Constant, p 3 , 
10‘ 3 C/(m 2 * C) 

— 

-0.25 


18 


Table n. - Convergence Studies of Beam Element 


Aspect Ratio 


[Og/p] Open Loop Case 
u w 


[0/90/±45] s 100 V Case 
u w 


400 

2.01475 e-4 

8.42061 e-3 

200 

2.01474 e-4 

8.42061 e-3 

100 

2.01471 e-4 

8.42061 e-3 

50 

2.01468 e-4 

8.42059 e-3 

25 

2.01468 e-4 

8.42059 e-3 


1.43711 e-4 
1.43683 e-4 
1.43645 e-4 
1.43632 e-4 
1.43631 e-4 


1.14907 e-2 
1.14906 e-2 

1.14905 e-2 

1.14905 e-2 

1.14905 e-2 
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Appendix A: Constitutive Equations: Constant 
Through-the-Thickness Displacement 


The simplified constitutive equations for constant through-the-thickness displacement are 


I s - 1 

► = 

0 

H 

1 



► + 

O 

J*- 
1 

I'M 

► + 

«1 

Ul 


O 

is* 


,°5 J 


[0 0 J 

1M 


0 


^3 “ ^31 a i + ®33 ^3 + Pi ® 

or in semi-inverted form 


(".1 

► = 

c n 

0 

P-i 

► — 

0 

*31 

1 

f E .] 

► — 

V 

\ J 


0 

^33 . 

I s * 1 


0 

0 


l*J 


0 


(A-l) 


(A-2) 


(A-3) 


J > 3 = e 31 5j + €33 E 3 + p 3 0 


(A-4) 
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Appendix B: Discrete Layer Laminate Matrices: 
Constant Through-the-Thickness Displacement 


The laminate density matrices are 

Pu = H / b P ♦*($ * 

/-I JZ 


*33 = E / ‘ b P 
1-1 Z 


the lamina te stiffness matrices are 


~ ]£ f b C 55 & 
m z 


4-E 

M J z * 


*>u f b C n V(z) I|r"(z) & 

1-1 Jz 

n fan . f h C J, 

D » * E J * C 55 -r 5 - * 


the piezoelectric laminate matrix is 

£,r - i: / » «,. ♦*« * 


(B-l) 

(B-2) 

(B-3) 

(B-4) 

(B-5) 

(B-6) 

(B-7) 
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the dielectric permittivity laminate matrix is 




dz dz 


the external forces are 

ft -it I b h **fe> * 

« z 

Si - E / b h & 

'*1 z 


4 


it 



J* $ i|r*(z) 


where L is the number of plies. 


(B-8) 

(B-9) 

(B-10) 

(B-ll) 
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Appendix C: Finite Element Matrices: 
Constant Through-the-Thickness Displacement 


The mass matrices are 



<> * J 

r P n R*(x) R J (x) dx 

(C-l) 

II 

r P 33 R‘(x) R j (x) dx 

(C-2) 

The stiffness matrices are 



K n = f ( D n dR ‘^ dRJ ty + D 55 R‘(x) R j (x) } dx 
Jx dx dx 

(C-3) 

- /, 

B k R‘(x) dx 

t dx 

(C-4) 

It 

E 31 /?'(*) dr 

i ctx 

(05) 

k »-L 

«<M «*'<*) * 
55 dx dx 

(C-6) 

Ml-- 

f x G 33 R*(x) R j (x) dx 

(07) 
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The external force vectors are 


f{ = U R'(X) ] x 

(C-8) 

Fj = [f 3 R*(x) } x 

(C-9) 

Q l = [q R*(x) ] x 

(C-10) 
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/ 1 Bottom Surface at 20°C 

Figure 1(b): [p/0/90/±45] s Beam Under A Thermal Gradient 
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Figure 2(a): [0,/p] Beam Mesh for Developed Finite Element Program 



Figure 2(b): [0,/p] Beam Mesh for Conventional Finite Element Program 
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Axial Displacement (m) 
(E-3) 



Beam Length (m) 


— Currant 

— FEM 


Figure 3: Axial Displacement of a [0,/p] Beam Under A Uniform 
Thermal Load (100 *C) in Closed Circuit Conditions 




Transverse Displacement (m) 
(E-2) 


0.80 
0.64 
0.48 
0.32 
0.16 
0.00 

0.00 0.05 0.10 0.16 0.21 0.26 
Beam Length (m) 

Figure 4: Transverse Displacement of a [0,/p] Beam Under A Uniform 
Thermal Load (100*C) in Closed Circuit Conditions 



28 




Axial Displacement (m) 
(E-3) 



Op«» Loop 
Closad Loop 


Beam Length (m) 


Figure 5: Axial Displacement Comparisons of a [QJp] Beam Under A 
Uniform Thermal Load (100* C) in Open and Closed Circuit Conditions 
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Transverse Displacement (m) 
(E-2) 



Open Loop 
Closed Loop 


Beam Length (m) 


Figure 6: Transverse Displacement Comparisons of a [0,/p] Beam Under A 
Uniform Thermal Load (100 *C) in Open and Closed Circuit Conditions 
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Electric Potential (Volte) 



Figure 7: Measured Sensory Electric Potential of a [Og/p] Beam 
Under A Uniform Thermal Load (100*C) 
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Transverse Displacement (m) 
<E-2) 



■■ Open Loop 

— - - OV 

— — — 100V 

• 200V 

300 V 

— 400V 

— — 500 V 


Beam Length (m) 


Figure 8:Active Compensation of the Transverse Displacement 
of a [<yp] Beam Under A Uniform Thermal Load (100 *0) 
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Axial Displacement (m) 
(E-3) 



Beam Length (m) 


Figure 9: Axial Displacement of a [p/0/90/±45], Beam Under 
A Thermal Gradient in Closed Circuit Conditions 
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Transverse Displacement (m) 
(E-1) 



0.00 0.05 0.10 0.16 0.21 0.26 

Beam Length (m) 

Figure 10: Transverse Displacement of a [p/0/90/445]. Beam Under 
A Thermal Gradient in Closed Circuit Conditions 
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Axial Displacement (m) 
(E-3) 



Figure 11: Active Compensation of the Axial Displacement of a 
[p/0/90/±45] g Beam Under A Thermal Gradient 


0 V 
100 V 
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0.00 

0.00 0.05 0.10 0.16 0.21 0.26 
Beam Length (m) 


Figure 12: Active Compensation of the Transverse Displacement of a 
[p/0/90/±45], Beam Under A Thermal Gradient 
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Electric Potential (Volts) 



Beam Length (m) 


Figure 13: Measured Sensory Electric Potential of a 
[p/0/90/±45]. Beam Under A Hernial Gradient 
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